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Abstract For a /3- Jacobi ensemble determined by parameters ai,a2 and n, under the restriction 
that the three parameters go to infinity with n and a\ being of small orders of ct2, we obtain both 
the bulk and the edge scaling limits. In particular, we derive the asymptotic distributions for the 



largest and the smallest eigenvalues, the Central Limit Theorems of the eigenvalues, and the limiting 
distributions of the empirical distributions of the eigenvalues. 

> 
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1 Introduction 

Let (3 > be a constant and n > 2 be an integer. A beta-Jacobi ensemble, also called in the 
literature as the beta-MANOVA ensemble, is a set of random variables (Ai, • • ■ , A„) S [0, 1]™ with 
^ r*| probability density function 

//W 3 (A) = ^" J] |A i -A/.;QAr P (l-A*) a2 - p , (1.1) 

l<i<j<n i=l 

where ai, a 2 > § (n — 1) are parameters, p = 1 + ^(n — 1), and 

S : cf ».« _ A r ( i + f)rK + a2 -f(,,-i)) 

The ensemble is associated with the multivariate analysis of variance (MANOVA). For (3 — 1,2 and 4, 
Q\ , the function //3(A) in (jl.lj) is the density function of the eigenvalues of Y* Y(Y*Y+Z*Z)~ 1 with ai = 

^nii and 0,2 — Qni2, where Y = Y mi xn and Z = Z m2 x „ are independent matrices with mi, 7712 > n, 



and the entries of both matrices are independent random variables with the standard real, complex 
or quaternion Gaussian distributions. See [14] and [57] for = 1,2. Other references about the 
connections between the Jacobi ensembles and statistics are [6| [12 | rT3] . fT4] . ri9| . f28l [34l [42 j [43 ^ [57]. 

In statistical mechanics, the model of the log gases can be characterized by the beta-Jacobi 
ensembles. A log gas is a system of charged particles on the real line which are subject to a 
logarithmic interaction potential and Brownian-like changes. If the particles are contained in the 
interval [0, 1] and are also subject to the external potential Y^i=i( I1 ^'~ 73) 1°§ Ai+( £ ^- — ^) log(l— Aj), 
where r = ^ai — n and s = ^0,2 — n, and (3 is the inverse of the temperature, then it is known that 
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the stationary distribution of the system of charges in the long term is the Jacobi ensemble as in 
(EH, see, e.g., [2] [23 E3 ES| . 

The beta- Jacobi ensembles also have connections to other subjects in mathematics and physics, 
for instance, lattice gas theory [32j[34], Selberg integrals [33l [55l [58] and Jack functions [4], 150], 162] . 

Now we briefly recall some research on the beta- Jacobi ensembles. Lippert [53] gives a model to 
generate the beta- Jacobi ensembles (see also [55] for a similar method used in the construction of the 
beta-circular ensembles). In studying the largest principal angles between random subspaces, Absil, 
Edelman and Koev [T] obtain a formula related to the Jacobi ensembles. Edelman and Sutton [59] 
study CS decomposition and singular values about these models. Dumitriu and Koev [22] derive 
the exact distributions of the largest eigenvalues for the ensembles. Jiang 02] derives the bulk and 
the edge scaling limits for the beta- Jacobi ensembles for (3 = 1,2 when p and a\ in (|1.1[) are of 
small orders of a 2 ■ Johnstone [49] obtains the asymptotic distribution of the largest eigenvalues for 
(3 = 1, 2 when oi, a 2 and p in (II. 1|) are proportional to each other. Recently, Demni [16] investigates 
the beta-Jacobi processes. 

In this paper, for the beta-Jacobi ensembles, we study the asymptotic distributions of the largest 
and smallest eigenvalues, the limiting empirical distributions of the eigenvalues, the law of large 
numbers and the central limit theorems for the eigenvalues. Before stating the main results, we need 
some notation. 

Let j3 > be a fixed constant, n > 2 be an integer, a\ and be positive variables. The following 
condition will be used later. 

nj3 

n — > oo, ai — > oo and — > oo such that a\ = 0(^/02), n = 0(^/02) and > 7 € (0,1]. (1-3) 



Then, for a sequence of probability measures {/i„; n — 0, 1, 2, • ■ ■ } defined on (R fe ,S(R fe )), we know 
fj, n converges weakly to (1q if and only if tZ(^t n ,/xo) — > as n — > oo, see, e.g., [21J. Similarly, we 
say that a sequence of random variables {Z n ; n > 1} taking values in R fc converges weakly (or in 
distribution) to a Borel probability measure \x on R fc if Ef(Z n ) — ► J Rk f(x) ji(dx) for any bounded 
and continuous function f(x) defined on M. k . This is also equivalent to that d(£(Z n ), /i) — * as 
n — > 00, where C(Z n ) is the probability distribution of Z n , see also [21] . 

For 7 S (0, 1], let 7 m i„ = (^7— l) 2 and 7 ma a; = (\/7 + I) 2 - The Marchenko-Pastur law is the 
probability distribution with density function 



2ai 



For two Borel probability measures fi and v on Mr, recall the metric 




(1.4) 





2tt~/x Vv^ Imin) (Ti 

0, 



i n ax 



x) , if x £ [7, 



otherwise. 



mm j 7maa:J ; 



(1.5) 
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The following is about the limiting distribution of the empirical eigenvalues of the beta-Jacobi 
ensembles. 



THEOREM 1 Let Ai, • • • , A„ be random variables with density function ff3 t ai,a 2 W as * n (GUP- Set 

1 - 

i=l 

for n > 2. Assuming \1.3\) , then d(/i n ,Mo) converges to zero in probability, where fio has density 
c ■ f-y(cx) with c = 2^/(3 and f 1 {x) is as in 1 1. 5\) . 

The next result gives the weak laws of large numbers of the largest and smallest eigenvalues for the 
beta-Jacobi ensembles. 

THEOREM 2 Let Ai, • • • , A„ be random variables with density function fp^ ai ,a 2 (-M as * n GUP- Set 
^max{n) = max{Ai, • • • , A„}, and A m i„(n) = min{Ai, • • • , A„}. Assuming U.3\) , we have that 

— ■ Amax{n) — > p- and — • X m in{n) -* p 

n 27 n 27 

in probability. 

Here is the central limit theorem for the eigenvalues for the model in 

THEOREM 3 Let Ai,--- , A n be random variables with density function ff3. ai .a 2 {^) as i- n {23P- 
Given integer k > 1, define 

' 1 1 fi\fi-l\ 

17' 



I (i\fi-l 
,r) \ r 



j=l r=0 

for i > 1, where c = 2^/(3 and 7 is as m il.3\) . Assuming il.3\) , then (X\, • • • , X^) converges weakly 
to a multivariate normal distribution -/Vfc(/z, S) for some \i and S given in Theorem 1.5 from 



Killip 51j obtains the central limit theorem for J27=i ^C^» e (°' wnere a < 6 are two constants. 
Theorem [3] is the central limit theorem for homogenous polynomials of Aj's. 

Thanks to the recent results of Ramirez and Rider [50] and Ramirez, Rider and Virag [55], we 
are able to investigate the asymptotic distributions of the smallest and largest eigenvalues for the 
beta-Jacobi ensembles next. Look at the operator 



Tp. a = - exp[(a + l)x + -j=b{x)] — |exp[-ax - —b(x)} 

where a > —1 and (3 > are constants, and b(x) is a standard Brownian motion on [0, 00). With 
probability one, when restricted to the positive half-line with Dirichlct conditions at the origin, 7g ja 
has discrete spectrum comprised of simple eigenvalues < Ao(/3, a) < Ao(/3, a) < ■ ■ ■ f 00 as stated 
in Theorem 1 from [60 . 

For a sequence of pairwise different numbers 01, • • • ,a n , let > > • • • > aS n ' be their 
order statistic. The following is the limiting distribution of the first k smallest eigenvalues. 
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THEOREM 4 LetXi,--- , A„ be random variables with density function fp ai a2 (X) as in $1.1]) . Let 
c > be a constant, and 2j3~ 1 a\ — n = c. If n — > oo and a 2 — > oo such that n — 0(^/02), then 
{2f3- 1 na 2 ) ■ {\ (n \- ■ ■ ,\ {n - k+1) ) converges weakly to (A (/3, c), A x (/3, c), • • • , A k -i(/3, c)). 

Now look at another random operator 

- H ' = £- X -& (L6) 

where b x is a standard Brownian motion on [0, +00). For AgR and function ip(x) defined on [0, +00) 
with ip{0) = and J °° ({ip') 2 + (1 + x)ifi 2 ^ dx < 00, we say (tp,X) is an eigenfunction/eigenvalue 
pair for — Tip if J °° %f) 2 (x) dx = 1 and 

^(i) = ^(i)ii + (i + A)^) 
holds in the following integration-by-parts sense, 



4>'{x) - V'(Q) = -JjjiK*)** + J --^b y i>'(y)dy + J (y + \)^(y)dy. 

Theorem 1.1 from [60] says that, with probability one, for each k > 1, the set of eigenvalues of — Hp 
has well-defined fc-largest eigenvalues Afc. Recall (|1.1|) . set 



(v^+^/S-^i and a n = - ^ H - LJ 



Our last result is about the limiting distribution of the first k largest eigenvalues of the beta-Jacobi 
ensembles. 

THEOREM 5 For each k > 1, let Afc &e t/ie /c-f/i largest eigenvalue of —Tip as in U.6]) . Let 
Ai, • • • , A n &e random variables with joint density function fj3, aua2 {X) as in 11. 1)) . Assuming il.Sty , 
then a n ((2a 2 ^ 1 )X W — uin),. , converges weakly to (Ai, • • • , Afc). 

REMARK 1.1 Dumitriu and Koev [22] derive the exact formulas for the cumulative distribution 
functions of the largest and smallest eigenvalue of the beta-Jacobi ensembles as in \1.1\) for fixed 
parameter j3, a,i,CL2 an d n - Here we complement their work by proving the asymptotic distributions 
in Theorems [7] and 

REMARK 1.2 In UW, Jiang study Theorems^ [H [3] and\3\for f3 = 1 and 2, which are special 
cases of the current theorems. The method used in JJ^j is the approximation of the entries of Haar- 
invariant orthogonal or unitary matrices by independent and identically distributed real or complex 
Gaussian random variables. 

Now, let us state the methodology used in our proofs. In fact, we employ a different approach 
than the standard ones in the random matrix theory. Some of the standard tools are the moment 
method in [7J [T51 (201 HH1 the Stieltjes transformations in [SH [55], the analysis to study the 
probability density functions of eigenvalues in [2 [46J [4JJ [65j [66j [67] , the large deviation method 
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for obtaining the law of large numbers in [5J [37] , and the application of the free probability 
theory in [TUl [TT] . A technique on refined estimates of the smallest eigenvalues is used in [M] . 

Another way to study random matrices is using the known conclusions, and connect them with the 
target ones by approximation. For example, large sample correlation matrices can be approximated 
by the Wishart matrices [45] ; a large dimensional, Haar- invariant matrix from the classical compact 
groups can be approximated by matrices with independent Gaussian entries |401 144j . The study of 
local statistics of the Wigner matrices with non-Gaussian entries can be approximated by Gaussian 
entries [55] . 

In this paper, we approximate the beta-Jacobi ensembles by the beta-Laguerre ensembles through 
measuring the variation distance between the eigenvalues in the two ensembles (Theorem[B]in Section 
[3J . Then the known results for the beta-Laguerre ensemble are used to get those for the beta-Jacobi 
ensembles. 

Now we would like to mention some future problems. Notice that all the theorems above are 
based on the restriction fL3j . We think it could be relaxed in some situations. One possible way is 
that, instead of using the uniform approximation in (|3.ip . one can treat case by case for the statistics 
concerned in the above theorems. For example, to improve Theorem [1] one could directly evaluate 
the moments Yl7=i for fc > 1 by computing the integration with respect to the density function 
fi3,a 1 ,a 2 (A) i n lll-ip . and then check what restrictions on a%, a 2 and n can make the integral close to 
the corresponding quantity in the Laguerre case. 

Two other problems are discussed in Remarks 13.11 and 13.21 

Finally we give the outline of this paper. In Section [2[ some known conclusions and some 
results on the beta-Laguerre ensembles are reviewed and proved, respectively. They will be used in 
the proofs for the beta-Jacobi ensembles. In Section [3J an approximation theorem for the Jacobi 
ensembles by the Laguerre ensembles is derived. In Section [4[ we prove the main results stated in 
this section. In Section [5j some known and useful results are collected for our proofs. 



2 Some Auxiliary Results on /3-Laguerre Ensembles 

Let P > be a constant, n > 2 be an integer, p = 1 + 4(n — 1) and parameter a > §(n — 1). A 
/3-Laguerre (Wishart) ensemble is a set of non-negative random variables (Ai,-- - , A n ) := A with 
probability density function 

n 

faaW=4' a II |A*-A/.n\ a_P -e-*^ A «, (2.1) 

l<i<j<n i— 1 

where 

" TV1 -u H\ 

^ = 2-™TT / b ■ ( 2 - 2 ) 

f = \T(l + lma-§(n-j)) 

One can see |25j for the construction of a matrix to generate eigenvalues with such a distribution. If 
X = (xij) is an m x n matrix with m > n, where Xij's are independent and identically distributed 
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random variables with the standard real normal (/3 = 1), complex normal (/3 = 2) or quaternion 
normal (/3 = 4) distribution, then /^(A) is the density function of the eigenvalues A = (Ai, • • • , A n ) 
of X*X with a = § m for /? = 1, 2, or 4. See [H [3j [57] for the cases (3 = 1 and 2, and [55] for 
/3 = 4, or (4.5) and (4.6) from [SI]. 

It is easy to see from Theorem A[S|in Appendix that the following is true. 

LEMMA 2.1 Let A = (Ai, • • ■ , A„) be random variables with the density function as in If 
n — > oo, a — > oo and n[3/(2a) — * 7 < 1, then 

i-l 



W ^ Z A J merges <0 (^) H ^ (*) (* r jY in probability; 
(**) ^ Z A J ~ (I) n Z ^ (*) (* r ^ ^ con« e r 5 e S to A a, 2 ) 



in distribution for any integer i > 1, where /i,; and a 2 are constants depending on 7, /3 and i onZy. 
A large deviation result in [23 [35] implies the following "law of large numbers" . 

LEMMA 2.2 Lei Ai, • ■ • , A n 6e random variables with the density as in \2.1\) . Assume n(3/(2a) — > 
7 € (0,1], and let 7 m j n = (^/7 — l) 2 and 7 maa; = (v/7 + l) 2 - Let /j, n be the empirical distri- 
bution of Yi := Ai7 / (n/3) for i = 1,2, ••• ,n. Then \x n converges weakly to the distribution floo 
with density f-y(x) as in il.5\) almost surely. Moreover, liminfn—too Y max (n) > jmax a.s. and 
limsup rl ^ 00 Y min (n) < j mm a.s., where Y max (n) = max{Yi, • • • , Y n } andY min (n) = min{Yi, • • • , Y n }. 

The first part of the above lemma is also obtained in [53] . 

Proof of Lemma 12.21 From (|2.ip . it is obvious that the joint density function of {Yi, 1 < i < n} 
is 

g( Vl , • • • , y n ) = Const ■ J] \y t - y/ f[ y*~* e~^^ v* 

l<z<j<n i — 1 

for yi > 0,- • • ,y n > 0. Take 7(71) = a - p, p(n) = n,Q(t) = f3t/(2j) in Theorem API] Then 
p(n)/n = 1 and j(n)/n = (a — p)/n /3(7 _1 — l)/2 as n — > 00 since p = 1 + (3(n — l)/2. According 
to Theorem Afl] /j, n satisfies the large deviation principle with speed {1/n 2 ; n > 1} and rate function 



l[v) = 

1 ' 2 



J f log \x - y\dv{x)dv(y) + / (^ - \ 1] logs) dv{x) + B 

= B + A ^-1. | y bg |z - y| di/(y) + ~ / (z - (1 - 7) logs) t&/(a:) 

where B is a finite constant. By Theorem Aj2[ the unique measure /Xqq on [0, 00) to achieve the 
minimum of I{v) over all probability measures on [0, 00) is the Marchenko-Pastur law with density 
function f 1 {x) as in (11. 5|) . 

For e > 0, let F = {v; d(u,Hoo) > e} f where d is as in (|1.4[) . Then, F is a closed set in the 
weak topology. By the large deviation upper bound, there exists a constant C > such that 
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P(d(fi n , ^oo) > e) < e~™ c as n is large enough. By the Borel-Cantelli lemma, d(fi n , [i^) —> a.s. 
as n — > oo. The first part of the conclusion is proved. 
For any integer fc > 1, it is easy to see 

/ y fc d/i„(y)) and > / d/i„(j/) ) . 

JO / imin\n) \Jo J 

Since /z„ converges weakly to almost surely, by the Fatou lemma, 

i/fc 



n— »oo 



for any integer fc > 1. Letting fc — > oo, we have 

liminf y moa; (n) > 7 ma x a.s. and lim sup Y min (n) < j min a.s. ■ 

n >oo n — 

In what follows, the notation x 2 ( s ) stands for the x 2 distribution with degrees of freedom s; x(s) 
denotes a positive random variables with (x(s)) 2 following the x 2 ( s ) distribution. 



LEMMA 2.3 Let X have the Gamma distribution with density function 

f{x\a,6) 



0, otherwise, 



where a > and > are constants. Given b > and > ; sei g(a) := P(X > 6). T/ien .g(a) is 
increasing over (0, oo). In particular, P(x 2 { s ) > b) < P(x 2 {t) > b) for any < s < t and b > 0. 

Proof. Since x 2 ( s ) nas probability density function f(x\a, 9) with a — p/2 and 9 = 2, we only need 
to prove the first part of the theorem. 

Obviously, X/9 has density function f(x\a, 1), without loss of the generality, we prove the con- 
clusion by assuming 9 = 1. First, noticing T(a) = J °° x a ~ 1 e~ x dx for a > 0, then 

f=o a — x J 

^ :=g(a+1)= /°°^e-x^ 
for a > —1. Use the fact that = x a logx for any x > and en e M to obtain that 

d<7i(a!) /fc 00 £ a (logx)e _a: dx ■ x a e~ x dx — x a e~ x dx ■ f^ x a (logx)e~ x dx 



da ( f °° i^e- 1 dx) 2 

-(/»(&) - h{0)), 



J °° i a e x dx- x a e x dx 



T(a + l) s 



where h(b) := / b °° x Q (logx)e x dx/ J b °° x a e x dx for 6 > 0. Fix a > — 1. To show 9l J^' > for all 
a > — 1, it is enough to show that /i(6) is strictly increasing on b e [0, oo). In fact, 
d/i(6) 



db 



(/~ x«e-- dxf 
> 



/OO /■ oo \ 

x a e- x dx + b a e- b J x a (\ogx) e - x dxj 
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since a;"(loga;)e x dx > (log 6) x a e x dx for all b > 0. ■ 

LEMMA 2.4 Let {a n ; n > 1} and /3 &e positive constants with 2a n > (n — l)/3 for all n > 1. 5et 
T„ = {(n - i + 1)/?, 2a„ - (i - l)/3; 1 < i < n} for n > 1. Le£ {x(i), f € T n ; rc > 1} &e a se£ o/ 
random variables defined on the same probability space. Define 

6 n = max I v(t) — ytl and p„ = max I Y 2 (t) — t| (2.3) 

teT„ teT„ 

/or n > 1. // n/3/(2a n ) — ► 7 € (0,1], </ien 5 n / y^n — * a.s. and Pn/n — * a.s. as n —» 00. 

Proof. First, noticing |T n | < 2n and max{t; t G T„} < nj3 + 2a n for all n > 1. Since n/3/ (2a„) — > 7 
as n — * 00, there exists a constant C > 1 such that max{i; t E T n } < Cn for all n > 1. Then, for 
any e > 0, 

P(<5„>Vne) < 2n ■ maxP (|x(i) - Vt| > V^e) 

< 2n • max P (\x(t) - Vi\ > y/ne] (2.4) 

0<t<Cn V / 

for all n > 1. Now, taking p n = [n 1 / 3 ] for n > 1, by Lemma [2~3l 

max P f |x(i) - V*| > V^e) < max P (x(t) > V^e - JpZ) 

0<t<p n \ / 0<t<p n 

< P (x(Pn) >Vn~t~ yfa) < P (\x(Pu) -VP^\> ^) 
as n is large enough since p n = o(y / n) as n — > 00. This and (|2.4p imply that, given e > 0, 

P (S n > V^e) < 2n ■ max P ( |x(i) — V*| > J (2.5) 

p-n<t<Cn V 2 y 

as n is sufficiently large. Now, the last probability is equal to 

p (x(*) > VI + ^) + ^ (x(*) < - ^) 

< p(x(W + 1)>^+^) +p(x(W)<^-^) (byLemmar 



< pf x(M + i)-y[^>-i + ^+^fx(W)-VW<i-^ 



< 2P max P(|v(fc) - Vfcl > — ) 
[t]<fc<[i]+i v ~ 3 

as n is sufficiently large, where the inequality \/T+ [i] — 1 < \/i < 1 + \/M f° r all t > 1 is used in 
the last step. Since |x(fe) - < \(x(k)) 2 ~ k\/Vk for any k > 1, by (|2~5| . 

P(d n >V^)<4». max pf l(x(*))'-fcl >^ (2 . 6) 
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as n is sufficiently large. Lemma 2.4 from [UJ says that 



P 



!ELife 2 -i)l 



for any n > 1 and c £ (0, y/n/2), where £i, £2, • ■ ■ , £n are i.i.d. random variables with £1 ~ AT(0, 1). 
Since > as n is large enough, the above inequality implies that 



max P 



iiicijjv _t 1 : — 

p„<fc<2C« V Vfc 



< max P 



^ l(x(fc)) 2 -fc| > VP~r 



liCd-A. J I , 

k<2Cn \ ^fk 



< 2e~ 



as n is sufficiently large. This and (|2.6|) tell us that, for any e > 0, P (<5„ > y 7 ™ 6 ) < e~™ 1/3 / 56 as n 
is large enough. Therefore, by the Borel-Cantelli lemma, 8 n /\/n — > a.s. asn-^ 00. 

Finally, since a 2 — o 2 = (a — b) 2 + 2b(a — b) for any a, 6 S K, recalling the constant C > 1 from 
(l2~4l). we have 



— < max h max 



«GT„ 



f2Vt \ x (t)-Vt\ 



< 



-p + (2C)--^^0 a.s 



asn^ 00. ■ 

By following the proof in [61j . we have a result below on the /3-Laguerre ensembles. It is also 
reported in Theorem 10.2.2 from [24] without proof. 

LEMMA 2.5 Let A = (Ai, • • • , A n ) be random variables with the density function as in t2.1]) . Set 
^max(n) = max{Ai, • • • , A„} and X m in = min{Ai, • • • , A„}. If n — > 00, a — > 00 and n(3/(2a) — > 7 S 
(0, 1], £/ien 

A m aa:( n ) o/i , / TT\\2 j -^mm(^) Q / 1 / TT\2 

-> p(l + V7 J a - s - arl « * P(l — VI J a - s - 



n 



asn-* 00. 



Proof. By Lemma \2. 21 it is enough to show 



Urn sup XmM <0(l+ Vt^ 1 ) 2 a-s. and 
hminf Amm(n) >(3(1- /F 1 ) 2 a.s. 



n— *oo 72 

From Theorem 3.4 in [25], we know the eigenvalues of P/3PJ have the same probability distribut 
as in (|2.1[) . where 



(2.7) 
(2.8) 



P« = 



X2a 

X/3(n-l) X2a~(i 



Xd X2a- P (n-1) ) nxn 
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where all of the 2n — 1 entries in the matrix are non-negative independent random variables with 
(Xs) 2 following the \ 2 distribution with degree of freedom s. It is easy to see that the first and last 
rows of BpBp are respectively are 

(XL, X2aX/3(n-l),0,- ■ -,0,0) and (0,0,0,- • • , XpXla-l3{n-2) , X/3 +xl a -0{ n -X))- 

For i = 2, • • • ,7i—l and n > 3, the «-th row is 

(0, ■ ■ • ,0, X/9(7i-t+l)X2a-(t-2)/3) Xp(n-i+l) + X2 a -(i-l)fi' X/J(n-i)X2a-(i-l)/9) 0) ' " ' > 0). 

s / 

zi/i position 

Recall the Gersgorin theorem: each eigenvalue of an n x n matrix ^4 = (ay) lies in at least one of 
the disks {2 £ C; |z — a^ l < \ a ij II f° r .7 = 1; 2, • • • , n. Then 

A ma x(n) < max{£/+, VF+; 2<i<n- 1} (2.9) 

where 

^ + = Xla + X2aX/3(n-l) ; ^ + = \% + X2o-/3(n-l) + X/3X2a-/3(n-2) and 
Wj = X/3( ra -i+l) + X2a-(i-l)0 + X/3(n-i+l)X2a-(i-2)/3 + X/3(n-i)X2o-(i-l)/3 ( 2 -10) 
for 2 < i < n — 1 , and 



where 



A min (n) > min{U-,V-,Wr; 2<i<n-l} (2.11) 



t 7 = X2a -X2aX/3(n-l), V = + X 2 a-/3(n- 1) ~ X/3X2a-/3(n-2) and 
W i = X/3(n-i+l) + X2a-(«-l)/3 ~ X/3(n-i+l)X2a-(i-2)/3 ~ X/3(n-i)X2a-(i-l)/3- 



Set 



Q+ i = /3(n-i + l) + 2a-(i-l)/3 + V/3(ti - i + l)(2a - (i - 2)/3) 



+ y/P(n-i)(2a-(i-l)(3) 
for 2 < i < n — 1 and n > 3. Reviewing the notation <5„ and T ra in Lemma \2A\ we have that 



|X/3(n-i+l)X2a-(i-2)/3 - VPip - i + l)(2a - (i - 2)/3) | 
= I (x/3(n-i+i) - \//3(n - i + 1)) [X2a-(i-2)() - y/2a-(i- 2)/3) 
+ v//5(n - i + 1) • (x2«-(i-2)/j - v/2a - (* - 2)/?) 

+ v/2a- (i-2)/3 (x/J(n-«-i) - >//3(n-i + l)) | 
< 5 2 n + (y<Jhl + V2^)8 n (2.12) 



uniformly for all 2 < i < n — 1 and n > 3. Moreover, 

lxjs(„_i+i) -/?("-< + 1)1 <Pn and |xL-(i-i)^-(2a-(i-l)/3)l < P« (2-13) 
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for all 1 < i < n and n > 1. Thus, 

I 9 max {W+} - max {Q+J| < max |W,+ - Q+J < tf(p„ + 5 2 n + 5 n ^) (2.14) 

2<t<n— 1 2<t<n— 1 ' 2<t<n— 1 ' 

as n is sufficiently large, if is a constant not depending on i or n. Evidently, < (in + 2a + 
2y/2afln — (y/pn + y/2a) 2 uniformly for all 2 < i < n — 1 and n > 3. By Lemma 12.41 and the 
condition that n/3/(2a) — > 7, we have 

max 2 < l <„-i{W / :+} /— Y2 
hmsup 5 — < p(l + V 7 J a - s - (2-15) 

n — >oo TL 

By the same but easier argument, the above is also true when max2<i< n -i{T4^ + } is replaced by U + 
and V + , respectively. Therefore, (|2.7|) is concluded from (|2.9p and 12. 151) . 

Now we prove ^Ml- Write (2a - (i - 1))/? = ~P + (2a - (i - 2))/3. It is easy to see 



Q~i : = /?(n - i + 1) + (2a - (i - l))/3 - y/P(n - 1 + l)(2a - (i - 2)/3) 

- v //5(n-z)(2a-( l -l)/3) 



> -/? + j9(n - i + 1) + (2a - (i - 2)/3) - 2^/5^ - i + l)(2o - (i - 2)/3) 
= -/3 + (V)9(n - i + 1) - ^2a - (i - 2)/3) 2 

for 2 < i < n — 1 and n > 3. Use the equality yfx — y/y = (x — y)j [yfx + y/y) to obtain 
{y/P{n-i + l) - y/2a-(i-2)fif = (P(n - 1) - 2af 

which is increasing in i G {2, 3, • • • , n — 1}. Since a > P(n — l)/2, we get 

(y/Kn-i + l)- V2a-^-2)/3) 2 = - ^ = *j = ^ 



mm 

2<i<n-l ' 



(V/9(n-l) + V^) 2 ' 

for n > 3. Combining all the above steps, and using the condition that nj3/(2a) — > 7, we arrive at 
i , mmf min 2 < I <„_ l{ Q- t} > /3(l-7-j 2 = _ ^ 2 ^ 
n^oo n (1 + v7~ 1 ) 2 

By (|2.12[) and (|2 . 13[) and the same argument as in (|2.14|) . we obtain 

I , min {W.r} - min {Q~ J| < max |WT - Q" .| < if (p n + 5^ + yft6 n ) 

2<i<n— 1 2<t<n — 1 2<z<n— 1 

as n is sufficiently large. This together with (|2.16[) and Lemma [2^41 concludes 
hminf ™2<<<n-i{WT} > 0(1 - yf^f a.s. 

n — >oo fi 

Similarly, 



U~ _ 2a- y/2a/3(n- 1) 



a.s. 



lim = lim = /?7 1 — fly/ 7 

n — >oo fi n — >oo fi 

V- 2a-/3(n-l) _! 
hm = lim = P7 — p a.s. 

n— >oo n n— »oo 77, 

Since 7 e (0, 1], it is obvious that min-^" 1 - Pypr 1 , fi^ 1 - p} > P(l - y/r^) 2 - By PTTTj) . the 
above three assertions imply (|2.i 



11 



3 An Approximation Result 



Let /x and v be probability measures on (M. m , B), where m > 1 and B is the Borel a- algebra on R™. 
The variation distance — v\\ is defined by 

\\n-v\\ =2- sup - = / |/(a;)-ff(aO|da;i---dsBn, (3.1) 
AeB ■/he™ 

if [x and have density functions f(x) and g(x) with respect to the Lesbegue measure. For a random 
variable Z, we use C(Z) to denote its probability distribution. The following is the key tool to prove 
the results stated in Introduction. 

THEOREM 6 Let /i = (/i l7 • • • , /i n ) and A ~ (Ax, • • • , A„) be random variables with density fp^ilA 
as in 12.1)) (taking a = a\) and fp,ai,aa{X) as * n (QI2P- Assuming 11.3)) , then \\£(2a2\) — £(/x)|| — > 0. 

REMARK 3.1 TTie condition lil.3)) is actually sharp for ||£(2a2A) — ^(a 1 )!! - * /or /3 = 1 and 2. 

In /act, Zei C/ be an N x N random unitary matrix with real ((3 = 1) and complex ((3 = 2), chosen 
with Haar measure. Decompose 

jj | A nit1l2 Cmx(JV-n2) \ 

\-S(W-ni) xn 2 -D(Af-ni)x(Af-n 2 )/ 

Assume rii > 712 a^rf n-x + n% < 2V. 5e£ y = A* ni n2 A ni ^ ri2 , we see from (3.15) in \3^j that the 
eigenvalues Ax, - • • , A„ 2 ofY have density function 

n 

Const- H \Xi-Xjf ■Y[Xt 1 - p (l-X i ) a2 - p , 

X<i<j<n i—1 

which belongs to the beta-Jacobi ensemble il.l)) with 

a% = 7^ n i, a 2 — 7j(N — ni), n — U2 and p = 1 + ~(n — 1). 

It is shown in \43tf for (3=1 and JJ^j f or (3 = 2 that the variation distance d! between yNA ni , n2 
and X n goes to zero if Tlx = o{yN) and n% = o(\f~N), where X n is the joint distribution of nxn2 
independent and identically distributed real Gaussian ((3 = 1) or complex Gaussian ((3 = 2) random 
variables. The orders n\ = o(yN) and 712 = o(y/N), which correspond to that ax = 0(^/02) and 
Tlx = o(y / a2 _ ) in kl.3)) . are also proved to be sharp for both cases in \42\ \4Stf - A further analysis shows 
that d — * if and only if ||£(2a2A) — £(/x)|| — ► 0. This tells us that the orders in Theorem^ are 
sharp for (3=1 and 2. However, it is not known whether the same remains true for other (3 > 0. 

REMARK 3.2 The condition "n(3/{2ai) — > 7" in \1.3)) is required in Theorem^ In the same 
contexts as in \43( for (3 = 1 and f42l for (3 = 2, the condition is not imposed. Although the 
condition is harmless in proving the main results stated in Introduction, it would be interesting to 
see if it can be removed for other (3 > 0. 

Now we start to prove the above theorem by developing several lemmas. 
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LEMMA 3.1 Let n > 2. Let fi — (pi, ■ ■ ■ ,(A n ) and A = (Ai,--- , A„) be random variables with 
density functions //3,oi(m) as * n CHIP (taking a = a\) and f/3 yai .a 2 W ^ n {ZHP > respectively. Then 

\\C{2a 2 \) - £( M )|| = £|tf„ • L„(/i) - 1| 

and E(K n ■ L„(/x)) = 1, where 

= a 2 M — —, ~R~\ fflW ( 3 - 2 ) 

£„(M) = e (1/2) E "=* " 4 • FT 1 - P-) ' J ( max Mi < 2a 2 ). (3.3) 

V Lao I l<i<n 

Proof. It is enough to show 

||£(2a 2 A)-£(//)|| = / \K n -L n {p)~l\ - fp, ai (ji) dp. (3.4) 

J[0,oa) n 

First, since p = %(n — 1) + 1, we have n(n — l)/5/2 + n(ai — p) + n = na\. It is easy to see that the 
density function of := 2a 2 A is 



0/3,01, a 2 (0) 



8,a\,a 2 



1 



n(ri-l)/3/2+ri(ai— p)+n n / \ a 2 — P 



l<i<j<n i=l V 



(, \ nai n / /j \ ai— p 

n nr- .-^ 
z/ l<i<j<n i=l v z ' 

for < i9i < 2a 2 and 1 < i < n, and is equal to zero, otherwise. Therefore, 



||£(2a 2 A) -£(m)II = / l3Aa!,a 2 (M) - //9,ai(i«)| rfy" 

J[0,oo)" 

,5/3,ai,a 2 (A*) 



[OjOo)* 1 //3,ai0) 



- 1| • fp, ai (^)dn. (3.5) 



Now, review //3,oi (a<) as in (|2.1() to see that 

m^uM = 2i ( j_\ nil- — i e st. m./2 

for < Hi < 2a 2 , i = 1, 2, • • • , n, and is zero, otherwise. It is easy to check that 

4' aua * / if 1 = 2 n ai f i \ M1 A r( ai + 0a -f(»-j)) 

c£' ai V2a 2 y l2a 2 ; 11 r(a 2 - f (n - j)) 

-nai TT r ( fl l + Q 2 - |i) 

Thus, 3/3,ai,a 2 (M)///3,ai (m) = ' £«(m)> which together with (|3.4p and (|3.5D yields the first conclu- 
sion. Finally, 

£(^„ • l„( m )) = | £ gf^ • Wm)^ = | gp, ai ,Mdix = i. ■ 
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LEMMA 3.2 Let h(x) = x logic for x > 0. For a fixed constant (3 > 0, an integer n > 1 and 
variables a\ > and 02 > 0, set b\ = -|ai and 62 = § a 2- -tf n — > 00, ai — > 00 and 02 — > 00 m a way 
ai = 0(^/02) and n = o( v /a2 _ ) 7 £/ien 

n / b — n\ 

+ 62 - * + 1) - fc(&a - * + 1)} = n&i ( 1 + log 6 2 + ) + o(l). 

Proof. Note that /i'(ir) = 1 + log a:, h"(x) — 1/x and /i^ 3 ^(x) = — 1/x 2 . Given x > 0, for any 
Am > —.To, by the Taylor expansion, 

h(x + Ax) - h(x ) = h'(x )Ax + l-h"(x a ){Ax) 2 + lh^(0(Ax) 3 

2 

= (l + logx )Ax+^-(AT) 2 - ^j(Ax) 3 
where £ is between xq and xo + Ax. Now take no = °2 — i + 1 and Ax = 61, we have that 

h(h + b 2 - i + 1) - M&2 - * + 1) 
= Ml + log(6 2 - i+ l)) + f^TT + o(|) (3.6) 
uniformly for all 1 < i < n. Obviously, 

log(6 2 - i + 1) = log 6 2 + log (\ - = log 6 2 - ^ + O (3.7) 

uniformly over all 1 < i < n as 

n — > 00. ai — > 00 and 02 — > 00 such that a± = 0(^/02), and n — 0(^/02). (3-8) 

Now, 

b\ 1 b\ b\ ( 1 1 



2 &a — i+ 1 26 2 2 V & 2-« + l 6 2 

_ _6f_ 6| i-1 
26 2 + 2 ' 6 2 (6 2 -i + l) 
6 2 /r;6 2 



26 2 V b 2 2 



uniformly for all 1 < i < n as (|3.8j) holds. Therefore, by (J3T6J) and (|3.7p . 

+ 62 — i + 1) — - i + 1) 
li l &i(*-l) 6 ? ^ A? + n 2 hi+n6? 

= 61 + 61 log 6 2 — 4 — - + -J- + o ' 



6 2 26 2 V bj 

uniformly for all 1 < i < n as (I3.8[) holds. Thus 



^ {&(&! + 6 2 - 1 + 1) - h(b 2 - i + 1)} 



6in(n - 1) n6j / 6f + n 2 6i + nb\ 



= nbi + nbi log 6 2 — h — + n ■ O , 

Zt>2 2t>2 \ 2 

, / , , bi-n\ bin / 6? + n 2 6 x + nb\ 

- w H 1+log62+ ^/ + ^ +n -°l^i — 1 

The conclusion follows because the last two terms are all of order o(l) as (13.8(1 holds. 
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K„ = exp<| (1 Q 7) f nS + o(l». ( :-!.!)) 



LEMMA 3.3 Let A'„ &e as m i fff.li)) . Assuming U.S]) , we then have 

-7)/ 
8a 2 7 

Proof. We claim that it suffices to prove 

^ = exp|^^M +0 (l)} (3.10) 

under assumption (|3.8p . If this is true, under the condition that nf3/(2a\) — > 7 < 1, it is easy to 
check that 



nai(ai - fn) _ (3 2 nH n (t n - 1) _ (1 - 7)/? 2 n 3 
2a2 802 8a27 2 



o(l) 



as (|3.8j) holds, where i n := 2ai/((3n) — > 7 x . Thus, (|3 - 9(1 is obtained 

f , 61 = !<zi and 6 2 = § 



Now we prove p.lOp . Set a = | , b\="^a\ and 62 = §02- It is easy to see that 



= " r(a(b 1+ 6 2 -i + l)) 

" * "to / f = \ r(a(6a-i + l)) ' 



Recall the Stirling formula: 



log r(z) = z log 2 — 2 — — log z + log \/2~k + — — h O 



as a; = Re (z) — > +00, where T(z) = J °° t z 1 e * d£ with Re (z) > 0, see, e.g., p. 368 from [35] or (37) 
on p.204 from [3]. It follows that 

logX„ (3.11) 
= — anbi log(a&2) 

+ ^ {a(6i + 6 2 - i + 1) loga(6i + b 2 - i + 1) - a(6 2 - i + 1) loga(6 2 - i + 1)} (3.12) 

- \ t log bl / +62 :; i +1 + O (3.13) 
2 ^ 6 2 - * + 1 \02-n) 

as (O holds. 

Now, write (an) log (ax) = (aloga)x + a(xloga;) and set h(x) = a; log a; for x > 0. Calculating 
the difference between the two terms in the sum of p,12j) . we know that the whole sum in (|3.12p is 
identical to 



a(log ajnbi + a ^(M & i + 62 - * + 1) - h(b 2 

i=l 

a(loga)n6i + anbi ^1 + log& 2 + + 

anbi + (anbi) log(a62) + anbi • ^ - + o(l) (3-14) 

269 
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by Lemma 15721 From the fact that log(l + x) < x for any x > 0, we have 

b\ + 62 — i + 1 v~s {., b\ \ nbi 



< 



b 2 - i + 1 r-f V 62 — + 1 7 62 — n 

2 — 1 Z — 1 

for any 62 > H. Thus, the sum of the three terms in (]3 . 13[) is equal to — omb\ + O (j^J a s (|3 
holds. Combining this and (|3.1ip - (|3.14p . we get 

logiT n = an&i • — h o(l) = * h o(l) 

2o 2 2a 2 

as {33]) holds. This gives l[3T0|) , ■ 

In the proofs next, we use op(l) to denote a random variable that goes zero in probability as taking 
a limit. 

LEMMA 3.4 Let \i = (pi, ■ ■ ■ ,fJ, n ) be a random variable with the density function as in \2.1\) with 
a = a\. Let L n {n) be as in h3.S\) . If il.S\) holds, then 

eXP \ 8a 2 y j"^)- 1 

m probability as n — > 00 . 
Proof. From (|3.3p , we see that 

n / \a 2 -p 

L„(M) = e (1/2) E ^ w • IT 1 - TT-) ' J ( max W < 2oa). 

■i— 1 x 7 

By Lemma O since ^ -> 7 € (0, 1] by (TO)) . 

maXl ^»^ -> 0(1 + v^) 2 (3.15) 
n 

in probability asn-> 00. Since n = 0(^/02), choose 5 n = (ny/a^) 1 ^ 2 , then S n /n — > 00 and <5 n /y / a 2 — > 
as taking the limit as in (|1.3p . Therefore, to prove the lemma, it is enough to show 

cxp( (1 ~ 7) f" 3 | ■L n (n) ^1 (3.16) 



8a 2 7 2 



in probability as n -» 00, where 



a 2 -p 



L^^ed/aJEr^Pi.JJfi-^.) . /( max Mi < «$„). (3.17) 

- L - L V 2d2 / l<i<n 

i— 1 x ' 

This is because, for any two sequences random variables {£„; n > 1} and {?7„; n > 1}, if — ► 1 in 
probability and P{£, n 7^ ?7n) -* as m 00, then rfo — ► 1 in probability asn-> 00. Rewrite 

n n ~j 

2 H Mi + ( fl 2 - p) 1 °g( 1 _ f 



1G 



on Q n :— {max!<^< n < S n }. Noticing log(l — x) = —x — (x 2 /2) + 0(x 3 ) as x — * 0, 

i=l z z i=l 2 i=l \ 2 i=l / 

on fl n . Now, on n n again, 

^Em?<^=(^) 3 -^--^0 (3.19) 

« 2 a 2 V V" 2 / V" 2 « 2 



as taking the limit in ([L~3|) . Recall p = 1 + §(n - 1). We have from (pTTg)) and ([3TTO|) that, on fi n , 

n 

(« 2 -p)X>(i-£) 

(a 2 - p)/3n 2 



«2 ~P | 




2a 2 1 




02 - P j 


( fi 2 n 3 


Sa 2 1 


{ 7 2 



2^ Mi 

i=l 



2a 2 7 

^ 2 ^3 



i=i 

as (jl.3[) holds. By Lemma T2.1[ as n — > oo, 



£f / 8a 2 7 2 V v/o2/ 



tE^ 4 -; -f>--n=>-iV(0,<7?); (3.20) 



n ~i 7 n t=i 7 

1 n fl 2 1 tf 2 1 ™ ft 2 

^^^(1 + -.2 7 )4(1| 7 ); - E A ~ ^(1 + 7)n => N(0, a 2 ) (3.21) 

i=l ' ' i=l ' 

P 

where a\ , cx 2 are constants depending on 7 only, the notation — * means "converges in probability 
to" and means "converges weakly to". Now, write (a 2 — p)/2a 2 = (1/2) — p/2a 2 , then 

a 2 - p ( I3n 2 -A \ 
fin 2 1 A pn if fin 2 A \ 

' i=l 

by (|3.20[) as (|1.3[) holds. Also, under the same condition, (a 2 — p)n 2 /a 2 = 0(n 2 /a 2 ) = o(l). It 
follows from (|3~2Tj) that 



»a 2 y 7 
(a 2 — p)7i 2 / fi 2 n 
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in probability as taking the limit in 1)1. 3|) . In summary, combining all the computations above, 

1 n n 
i—l i—l 

/3n 2 {a 2 -p)(3n 2 /3 2 n 3 (a 2 -p) 

= 17"^^ 8^1^ (1 + 7) + op(1) 

Ppn 2 (3 2 n 3 (3 2 n 3 p 
= 2^7"8^ (1 + 7) + W ( +7) + ° P( } 
on O n . Now, since p — 1 + |(ra — 1), n/^faa — > 0, we see that 

f3pn 2 I3 2 n 3 (3 2 n 3 p 

o = 1 h o(l) and 2 (1 +7) -> 

2a 2 7 4a 2 7 8a^7 2 

as (11.3[) holds. Thus, on f2„, 

n n 

5 X> + <« 2 -P)I>g<i-^) 

2—1 i—l 

P 2 n 3 P 2 n 3 | | ^ _ (7- l)P 2 n 3 | ^ 
4a27 8ci27 2 8a27 2 

as taking the limit in (| 1 .3ft . By reviewing (|3 . 1 7[) . we conclude (|3.16p . ■ 

The following is a variant of the Schcffc Lemma . 

LEMMA 3.5 Let {X n ; n > 1} be a sequence of non-negative random variables. If X n — ► 1 in 
probability and EX n — > 1 as n — > 00, then E\X n — 1 1 — s- osn-> 00. 

Proof. By the Skorohod representation theorem (see, e.g., p. 85 from [H]), w.l.o.g., assume X n — ► 1 
almost surely asrn 00. Thus, for any K > 1, we have X n I(X n < K) — ► 1 almost surely as?n 00. 
This gives that E\X n I{X n < K) — 1| — > and EX n I{X n < K) — > 1 as n — > 00. The second assertion 
and the condition that EX n — > 1 imply EX n I(X n > K) — > as n — > 00. Therefore, 

£|X„ - 1| < £|X n J(X„ < A') - 1| + EX n I(X n >K)^0 

as n — ► 00. ■ 



Proof of Theorem [6] It is known from Lemma T3. II that 

||£(2a 2 A) - £( M )|| = £|A„ • L n (p) - 1|. 
with E(K n ■ L n (fi)) = 1 for all n > 2, where /i has density fp,ai(p>) as m (BUI)- By Lemmas 13.31 and 

El 

A» = exp| v + o(l) | and exp | J-L»(p)-*1 

in probability as taking the limit in (|1.3|) . These imply that K n ■ L n (pL) — > 1 in probability as taking 
the same limit. Then the desired conclusion follows from Lemma l3~5l ■ 
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4 The Proofs of Main Results 



By using Theorem [6] developed in Section 02 we now are ready to prove the results stated in Intro- 
duction. 

Proof of Theorem [Q Set 



1 " 



n ■ i 



for (X[, • • • , X' n ) £ [0, +oo) n . Then, recall the definition of d in (|1.4|) . by the triangle inequality, 
|d(iU n ,/xo) - d(y n ,fio)\ < d{fi 

1 - 

= sup \~V] (f(n~ 1 a 2 X i ) - f{n' 1 a 2 X' i )) \ 
\\f\\ B L<i n t^i 

d 2 . , . 

< — • max Aj — A, , 

Tl l<i<n 

where the Lipschitz inequality \f(x) — f{y)\ < \x — y\ is used in the last step. This says that 
d(fx n ,fio)t as a function of (Ai, • • • A n ), is continuous for each n > 2. Thus, for any e > 0, there exists 
a (non-random) Borel set Act™ such that {d(fi n , fio) > e} = {(Ai, • • • , A„) £ A}. Then, by the 
definition of the variation norm in (|3.ip we see that 

P(d(M«, Mo) >e)< P{d(ji! n , fi ) >e) + \\C(2a 2 X) - C{fi)\\ (4.1) 

for any e > 0, where fi' n = (l/n)5Z"=i -W(2n) and fi = Qui, • • • ,/i„) has density fp, ai (p) as in (|2~T) 
with a = a\ and n/3/2ai — > 7 G (0, 1]. By Lemma 12.21 with probability one, 

1 " 

— > JfiT converges weakly to (4.2) 

n z — ' ""3" 

with density f^(x) as in Crg) . Write M,/(2n) = (^7 /n^cT 1 , where c = 2 7 //3. Then, by P~2"|) . 
with probability one, /x^ converges weakly to /io, where ^0 has density function c • f 1 (cx). Equiva- 
lently, d(fi' n , fio) — > almost surely. This, (14. ip and Theorem[S]in Section [3] prove the theorem. ■ 

Proof of Theorem [2j First, A max (n) and A m j„(n) are continuous functions of Ai, • • • , A n for any 
n > 1. Then 

p (,a 2 /3(1 + V7) 2 , . 

P — X m ax(n) — * > e 

\ n 27 

< P f 1 j-/w(n) - ^ t ^ 1 > e ) + ||£(2oaA) - £(fi)\\ (4.3) 



2n 27 

for any e > 0, where fi = (fi\, • • • ,/x n ) has density fp^if 1 ) a s in (|2.ip with a = a% and nj3/2a\ — > 
7 G (0,1]. From Lemma we know /i roax (n)/(2n) -> + Vt^ 1 ) 2 / 2 = £(1 + V7) 2 /( 2 7) in 
probability. This together with (|4.3[> and Theorem [6] in Section [3] yields the desired conclusion. By 
the same argument, (a 2 /n)X m i n (n) converges to (3(1 — \/l) 2 / (2l) in probability. ■ 
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Proof of Theorem [3] Evidently, 

\P((X U - ■ ■ ,X k ) e A) - P((Y U - ■ ■ 7 Y k ) e A)\ < \\£(2a 2 X) - 
for any Borel set iel', where 

C \ l x-^ 1 ( l\ fl—1 



j=l r=0 

7 \ ^-^ 1 / 1 \ ft — 1 



7 



for i > 1 (since c = 2"f/f3), and /i = (/ii, • • • , fJ- n ) has density fp^iilA as in (|2.1[) with a = a\ and 
n(3/2ai — > 7 S (0, 1]. The conclusion then follows from this, Theorem AJ3]and Theorem [5] in Section 
1 ■ 

Proof of Theorem[4l The assumption that 2j3~ 1 ci\ — n = c and n = 0(^/02) imply that n(3/2a\ — > 1 
and ai = 0(^/02). Thus Theorem [()] in Section [3] holds. 

Let • ■ • , #n) have density //5. ai as in (|2.ip with a = ai. Noticing, "/3Ai" and "a" in Theorem 
1 from [60j correspond to and "c" here, respectively. By Theorem 1 from |60) , for fixed integer 
k > 1, 



/3 ' 



,^6» (n fe+1) ^ converges weakly to (A (/3, c), • • • , A k -i((3, c)) 

as n — ► 00. By Theorem [5] in Section [31 

P((2a 2 A 1 ,--- ,2a 2 A„) € £„) - P((0 X , • • • ,0„) € £„) -> 

for any Borel set B n c K" for n > 1. From the Weyl perturbation theorem, we know that \^ is a 
continuous function of (Ai, • • • , A„) for any 1 < i < n. Combining the above two limits we obtain 

^ A («),..^ A M+i)j converges weakly to (A (/3, c), • • • , A k _i(/?, c)) 
as n — > 00 and a 2 — > 00 with n = c^-^/a^)- The proof is complete. I 
Proof of Theorem [5j Recalling ([HE]), lct 

Let (6*i, •• • , 6*„) have density //3 j(1 as in (|2.1[) . Noticing, "/3A 2 " in Theorem 1.4 from [515] corresponds 
to here; "k" in Theorem 1.4 from [59] is equal to 2/3 _1 a, and <|(n — 1) + 1 = p, and the fc-th 
lowest eigenvalue of Tip is the (n — k + l)-th largest eigenvalue of —Tip. Then by Theorem 1.4 from 
[59], 

'Q{1) _ 



converges weakly to (Ai,---,Afc) (4.5) 

!=!,— ,k 
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as n — > oo and a — ► oo such that n/a converges to a nonzero, finite constant. In other words, 

for any Borel set icM'. By Theorem [6] in Section [3l assuming (| 1 . 3[) and a = ai, 

P((2a2Ai,--. ,2a 2 X n )eB n )-P((9 u --- ,8 n )eB n )^0 

for any Borel set B n C M n for n > 1. The Weyl perturbation theorem says that <?(a;) := x"\ the 
i-th largest one in {x\, • • • , cc m }, is a continuous function of (xi 7 • • • , x m ) € R m for any 1 < I < m. 
Replacing a by ai in (|4.4p . the above two assertions conclude that 



er„ [ m„ converges weakly to (Ai, • • • , AfeJ. 

i=i.- ,k 



5 Appendix 

Let Q{x) be a real continuous function on [0, oo) such that for any e > 

lim xe~ £Q( - x) = 0. (5.1) 

X — > + oc 

For each integer n > 1, let p(n) be an integer depending on n. Let Ai, • • • , X p r n -\ be non-negative 
random variables with joint probability density 

/ P(") \ P(") 

exp -n^g(A,) n A " H II ^ X ^ ^ 

™ \ z=l y i=l l<i<j<p(n) 

where (3 > is fixed and 7(71) > depends on n. Let /x„ be the empirical probability measure of 
Ai, • • • , Ap(„). Under the weak topology, the large deviations for {/i„} is given below. For a reference 
of general large deviations, one can see, e.g., [TS] and [T7] . 

THEOREM A.l (Theorem 1 in JB3$ ) Assume p(n) / 'n —> re e (0, 1] and j(n)/n — > r > as n — > oo. 
TTien i/ie finite limit B :— lim,^,^ n~ 2 log Z„ exists and n > 1} satisfies the large deviation 

principle with speed {n~ 2 ; n > 1} and good raie function 

I(v):=-^y J J \og\x-y\dv{x)dv(y) + k, J (Q(x) - rlogx) du(x) + B (5.3) 

/or aZZ probability measure v defined on [0, oo). Moreover, there exists a unique probability measure 
vq on [0,oo) such that I(vo) = 0. 

In Theorem 1 from [38] or Theorem 5.5.1 from [37], the limit r above is required to be strictly 
positive. However, after a check, it is found that the conclusion also holds for r = 0. 
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THEOREM A. 2 (part of Theorem 8 from [38]) Let 7 € (0,1], j min = (1 - ^/l) 2 and ^ max = 
(1 + ^/7) 2 • For probability measure v, define 

J(v) = ~y J J lag\x-y\dv{x)dv(y) + '~ J (x - (1 - 7) loga;) dv{x). 

Then the unique minimizer of J{y) over all probability measures on [0, +00) is the Marchenko-Pastur 
law with density function / 7 (x) as in \1.5\) . 

Proof. Take a in [38J equal to 7. Notice 



V47 - (x - 1 - 7) 2 = yj(x - (1 - V7) 2 ))((l + ^) 2 - x) = yfe - T^Kw - ^) 
for all x G [7mm, 7moi]- Also, over all probability measure v on [0, 00), taking the minimum for 

I{v): = J J \og\x-y\dv{x)dv{y) + ^ J {x - {1 - i)\ogx) dv{x) 

-\ ( 3 7 ~ 7 2 log 7 + (1 - if log(l - 7)) 

is the same as doing so for J{v), where OlogO := as the convention. Then the conclusion follows 
from Theorem 8 in 1381. ■ 



THEOREM A. 3 (Theorem 1.5 from ]23}j ) Let Ai,-- - ,A n be random variables with the density 
function as in &2.1\) . Assume nf3/(2a) — > 7 < 1. Define 



*-S(^)'-g^C)C;> 



/or i > 1. Then, as n — > 00, (.Xi, • ■ ■ ,-X)fe) converges weakly to a normal distribution Nk{fx, S) /or 
some /1 and S. 

Acknowledgement The author thanks Professor Ioana Dumitriu for helpful discussions on the 
beta ensembles. 
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